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NONPERTURBATIVE RELATIVISTIC APPROACH TO PION FORM FACTOR: 
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Some predictions concerning possible results of the future JLab experiments on the pion form 
factor F^(Q^) are made. The calculations exploit the method proposed previously by the authors 
and based on the instant-form Poincare invariant approach to pion considered as a quark-antiquark 
system. Long ago, this model has predicted with surprising accuracy the values of Ftt{Q^) measured 
later in JLab experiment. The results are almost independent from the form of wave function. The 
pion mean square radius (r^) and the decay constant also agree with experimental values. The 
model gives power-like asymptotic behavior of F„(Q^) at high momentum transfer in agreement 
with QCD predictions. 
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I. INTRODUCTION 



The recent high accuracy experiments on the measure- 



ment of the pion form factor in the range of Q up to 
2.45 GeV2 P, i] (Q^ ^ -q^^ q- momentum transfer) 
and future JLab experiments up to « 6 GeV^ 01 
enhanced the interest to theoretical description of pion 
at high . 



It is usually believed that these future experiments 
will provide a meaningful test of the transition between 
perturbative and non-perturbative regions which is ex- 
pected at much lower in the case of pion that of other 
hadrons, in particular, of nucleon. At the present time 
different theoretical approaches to the pion form factor 
Ftt{Q^) exist. They are partly listed and described in 
Ref. [2] (section IV) (see also jSj). In the frameworks of 
some of these models, a certain agreement with existing 
experimental data is obtained for soft F-^(Q^). As to the 
region of high momentum transfer, the theoretical results 
differ from one another to a great extent. It seems us 
that the situation is such that one has almost no hope to 
find the appearance of perturbative degrees of freedom in 
the future JLab experiments on F^{Q'^). It is difficult to 
imagine that in the wide band of non-perturbative the- 
oretical curves there would not occur any one agreeing 
with the experimental data. Our opinion is that the large 
variety of non-perturbative predictions for future data for 
-Fir(Q^) makes it necessary to formulate the problem of 
detecting of perturbative effects in a slightly different way 
than it is usually done. Namely, we propose to accept 
one of the theories which describes correctly the exist- 
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ing data and continue the calculations for higher Q^. If 
it will occur that beginning from some values one needs 
to adjust the calculations by introducing the quark mass 
dependence on to agree the future data, then we will 
identify these values with the appearance of perturbative 
effects. It is natural that in the present paper we choose 
our own approach Q as an example for the demonstra- 
tion of the proposed scenario. 

The reasons for this choice are the following. The main 
one is that our approach has already demonstrated its 
predictive power: without any additional tuning of pa- 
rameters, we predicted in ReL M the values of F.,r{Q'^) 
obtained later in experiment 0, HHI- At the same time, 
the approach gives the correct values of the mean square 
radius (MSR), the decay constant /^r and the power- like 
asymptotic behavior. Certainly, other criteria of discrim- 
ination the approach may exist. For example, one can 
consider as a "correct" one an approach which gives a 
consistent treatment of the pion form factor in space-like 
and time-like regions. 

In the present paper, we use the approach to the pion 
form factor -FVlQ^) proposed in our papers about ten 
years ago [S] . Our approach presents one of the versions 
of the constituent quark model (CQM). The method is 
based on the dispersion approach to the instant form of 
the Poincare invariant quantum mechanics [lol | (see also 
the detailed version ll] and the review )■ 

Based on this approach and on the experimental data 
on the measurement of F^((5^) in the range of up 
to 0.26 (GeV)2 [H, we have obtained in 1998 the model 
function for the pion form factor for the extended range of 
higher momentum transfers [Q]. The experimental data 
obtained later 0, H, Q (see also the review of all experi- 
mental results in Ref. [2l] and references therein) for the 
range of larger by an order of magnitude coincide pre- 
cisely with our theoretical curve of 1998 [1] with no need 
of any additional fitting. This means that it is possible 
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to consider our calculations [6| as an accurate prediction 
of the present experimental data for the pion form factor. 
The model describes correctly the pion MSR and the de- 
cay constant too. It is important to notice that the 
dependence of our results for Ft^{Q^) on the form of wave 
functions is very weak Q. Moreover, our approach gives 
the correct power-like asymptotic behavior of Ft^{Q^) at 
oo. So, the model works well at high as well as at 
low values of Q^. 

Taking into account these advantages of our approach, 
it seems natural to hope that the model will continue to 
give a good description of experimental data at higher 
momentum transfers, in particular of the future JLab 
measurements in the range 2.45 < < 6 (GeV)^ (af- 
ter having withstand the test of tenfold increasing of 
range it may withstand another much smaller increase). 
If it will occur that the experimental data does not follow 
our theoretical curve then we shall adjust the theory by 
taking into account the quark-mass dependence on Q^. 
Within our approach, this depedence is a manifestation 
of appearance of the perturbative degrees of freedom. 

The paper is organized in the following way. We start 
in Sec. II with a brief review of the basic theoretical for- 
malism of our approach. The results of calculations and 
the comparison with the experimental data and other 
theoretical models are given in Sec. III. In Sec. IV the 
asymptotic behavior of the form factor is considered. Fi- 
nally, our conclusions are given in Sec. V. 



II. THE MODEL 

Our method is a version of the instant form of the 
Poincare invariant constituent-quark model (PICQM), 
formulated on the base of a dispersion approach (see, 
e.g., d, As is well known the dispersion approach 

is based on the general properties of space and time and, 
therefore, is to a certain extent "model independent". 
That is why the calculation of electromagnetic form fac- 
tors using dispersion approach are of distinguished char- 
acter as compared with other approaches. This advan- 
tage of our method is emphasized in Ref. [MJ (see, how- 
ever, the footnote [51||). 

The main point of our approch is the construction of 
the operator of electromagnetic current which preserves 
Lorentz covariance and conservation laws in the rela- 
tivistic invariant impulse approximation (so called mod- 
ified impulse approximation (MIA)) (lo| . This approxi- 
mation is constructed using dispersion-relation integrals 
over composite-particle mass, that is over the Mandel- 
stam variables s,s' [ll|. This variant of dispersion ap- 
proach was developed in Refs. [H, M, M, [S [H, S El 
and was fruitfully used to investigate the structure of 
composite systems. 

Let us recall some principal points of our approach 
[li I'^ variant of PICQM, pion electromagnetic 



form factor in MIA has the form 

FAQ^) = I d^/id\/7^(fc)go(s,g^s')¥'(^') • (1) 



Here tf{k) is pion wave function in the sense of PICQM, 
go{s,Q^,s') is the free two-particle form factor. It may 
be obtained explicitly by the methods of relativistic kine- 
matics and is a relativistic invariant function. 

The wave function in ([T]) has the following structure: 

Here M is the mass of the constituent quark. Below for 
the function u{k) we use some phenomenological wave 
functions. 

The function gQ{s,Q^ , s') is written in terms of the 
quark electromagnetic form factors in the form 



5o(s,Q^,s') 



(s + s' + Q2)Q2 



[X{s,-Q\s')f'^ Vl + 0V4M2 



x{{s + s' + Q^)[G%{Q^) + G%{Q^)] 



X cos + — e(5, Q\s'){GUQ^) 



+G«,(Q2))sin(c^i+^2)} . (2) 
Here A(a, 6, c) = + 6^ + - 2{ah + ac + be), 

ui and LU2 are the Wigner rotation parameters: 

c(s,Q^s') 



LUi = arctan 



UJ2 = arctan ■ 



A/[(Vi+V7)2+Q2]+^(yj+y^)' 



a{s,s')as,Q^s') 



M{s + s' + Q^)a{s,s') + Vss'{4:M^ + Q^)' 



a{s, s') = 2Af ^ + V7, 9{s, Q^, s') = d{s' - si) 
i?(s' — S2), i? is the step function, 



si,2 = 2Af2 + ^(2M' + Q')(s - 2Af2) 



T^VQHQ' + 4M2)s(s - 4AP). 

Note that the magnetic form factor contribution to 
Eq. ^ is due to the spin rotation effect only [23 |. 
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Here G^\[{Q'^) are electric and magnetic form factors 
of quarks, respectively. 

Let us note that we introduce electromagnetic quark 
form factors, in particulary, in order to obtain a descrip- 
tion of the maximal set of experimental data on pion, 
including the MSR and the decay constant simultane- 
ously, at the same values of the parameters of the model 

m. 

We use the following explicit form of the quark form 
factors: 



NpL = 16 



u{r) = Nt exp(-ar^/^ - (ir) 



a=ly/2M^, (3 = Mrb, Nt 



(5) 



where Cq are quarks charges and Hq — anomalous mag- 
netic moments which enter our equations through the 

sum Sn = Ku + K J. 



l + ln(l + (r2)02/6) 



(3) 



where (r^) is the quark MSR. 

Let us discuss in brief the motivation for choosing the 
explicit form ([3]). One of the features of our approach 
is the fact that the form factor asymptotic behavior at 
^ 00, M ^ does not depend on the choice of the 
wave function in Eq. ([T]) and is defined by the relativistic 
kinematics of two-quark system only I6l|. In the point- 
like quark approximation (Kg=0, {r'^)= 0) the asymp- 
totics coincides with that described by quark counting 
laws [13, (see also the recent discussion in 26]): 
-p7r(Q^) ~ ■ The asymptotic behavior of the pion 
form factor was considered in ^271,1281 (see also Ref. [29|]). 
The form ([3]) gives logarithmic corrections to the power- 
law asymptotics , obtained in QCD. So, in our approach 
the form ([3]) for the quark form factor gives the same 
asymptotics as is predicted by QCD. Let us note that 
another choice of the form of q uark form factor, for ex- 
ample, the monopole form [30|, changes essentially the 
pion form factor asymptotics so that it does not corre- 
spond to QCD asymptotics anymore. 

To calculate the pion form factor we use wave func- 
tions of different forms: harmonic oscillator wave func- 
tions (analogous to those used in the seminal paper [3l| 
continued recently in ^3^]), power-law- type wave func- 
tions with the explicit form motivated by perturbative 
QCD calculations at high [H, [s^, H^, wave functions 
with linear confinement and Coulomb-like behavior at 
small distances [3^. These functions are of the form: 

u{k) ^ NHoe^p{~k^/2b^) , 



5 2 



N{a,P) =9aV2ar - li^i -, ^ 



6' 3' 



-2^/3r(- ,FJ-,-,t 



7 4 



,345 
6/3^2^2 ( 1, 2' 3' 3' t 



t 



27 a2 



(6) 



where a , b are the parameters of linear and Coulomb 
parts of potential, respectively, Mr is the reduced mass 
of the two-particle system, b = (4/3) as, as =0.59 on the 
scale of the light mesons mass, pFq are hypergeometric 
functions, r(a;) is the Euler gamma function. 

The parameters of the model are the same as in Ref. @ 
where the motivation of the choice is described in detail. 
Let us note that for the constituent-quark mass M = 0.22 
GeV the values of parameters ([4]) - ([6l) were chosen in 
such a way as to ensure the pion MSR within experimen- 
tal uncertainties {rl)lLl = 0.657±0.012 fm ^ as well 
as the best description of the decay constant f-^exp = 
0.1317± 0.0002 GeV [S^- The sum of the quark anoma- 
lous magnetic moments is Sq = 0.0268, the quark MSR is 
(r^) ~ 0.3/Af2. The values of other parameters are 
the following: in the model (H b =0.3500 GeV (the 
decay constant is /^r = 127.4 MeV); in the model ([5]) 
n = 3, 5 =0.6131 GeV (/^ = 131.7 MeV); in the model 
® a =0.1331 GeV2 (/^ = 131.7 MeV). An interesting 
feature of our results is the fact that at the fixed con- 
stituent quark mass, the dependence of the pion form 
factor on the choice of the model ([1]) - (O is rather weak. 
The curves calculated with different wave functions but 
one and the same quark mass form groups 0]. From 
the theoretical point of view this weak dependence of 
our calculations on the model is the consequence of the 
dispersion-relation base of the approach. 



III. RESULTS OF THE CALCULATIONS 



Nho - 



(4) 



u{k)=NpL{k'lb' + n = 3. 



The results of the calculation of the charge pion form 
factor using the wave functions (|1]) , ([5]) , ([5]) and the value 
of constituent-quark mass M = 0.22 GeV (this parameter 
has been fixed as early as in 1998 @ from the data at 
Q2 < 0.26 (GeV)2 [l||) are shown in Figures [T] and H 
Let us note that our relativistic CQM describes well the 
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FIG. 1: (color online). Our predictions for the pion form fac- 
tor given in 1998 (full red lines; upper line: wave functions 
Q, lower line: wave functions ([SJ with n — 3 and ^) com- 
pared with data and with some other models. Green crosses 
represent data points of Amendolia et al. l^M- Other data 
points (all taken from Ref. 0|) are: reanalyzed points of Ack- 
erman et al. [s^ (full circle); reanalyzed points of Brauel 
et al. [3^ (full triangle); JLab results (full diamond, blue 
and yellow squares). Other theoretical curves are: QCD ap- 
proximation of Maris and Tandy !40] (dotted); perturbative 
QCD (leading and next-to-leading order) by Bakulev et al. 
[SI (green dash-dot-dotted); Nesterenko and Radyushkin [i^ 
(magenta dash-dotted); dispersion approach of Donoghue and 
Na [431 (blue dashed); holographic approach of Grigoryan and 
Radyushkin [13] (thick grey). 



experimental data for the pion form factor including the 
recent points ^ . The upper of our curves corresponds to 
the model (jl]), the lower - to the models ([5]) with n — 3 
and ([U, which lie close to one another. 

Let us emphasize that the parameters used in our cal- 
culations were obtained from the fitting to the experi- 
mental data up to ~ o.26 GeV^ [l3|. At that time the 
data for higher was not correlated in different exper- 
iments and had significant uncertainties. The later data 
for pion form factor in JLab experiments up to =2.45 
GeV^ were obtained with rather good accuracy. All ex- 
perimental points obtained in JLab up to now agree very 
well with our prediction of 1998. 

Let us discuss briefly the basic moments that pro- 
vide good results for the pion form factor in our ap- 
proach. First, throughout the calculation the condition 
of Lorentz-covariance and the conservation laws for the 
operator of electromagnetic current were satisfied. Sec- 
ond, the condition of accurate description of the pion 
MSR constrains the behavior of the wave functions in 
momentum representation (l4|), ([5]) at small relative mo- 
mentum of quarks or of the wave function in coordinate 
representation ^ at large distances because of special 



FIG. 2: (color online). Comparison of our predictions with 
other CQMs. Data and our curves are the same as in Fig. [1] 
Other theoretical curves are: those by C.-W. Hwang [i^] 
(blue dashed); Cardarelli et al. and (precisely coinciding 
with it) instant-form predictions 46] (magenta dash-dotted); 
Ref. |47| (green dotted). Predictions of an upgraded version 
of a seminal paper [sj] (Ref. [32] . Fig. 5, iriq = 0.22) coincide 
precisely with our upper curve. 



properties of the integral representation ([T]). Third, the 
condition of the best description of the decay constant 
defines constraints for the wave function at large relative 
momenta because the contribution of small relative mo- 
menta to the decay constant is suppressed as can be seen 
from the relativistic formula (see, e.g. [1,13): 



MV3 



fc2 dk 



(fc2 + M2)3/' 



■uik) 



So, our way of fixing the model parameters constrains 
effectively the behavior of wave functions both at small 
and at large relative momenta. The structure of our rel- 
ativistic integral representation ([T]) is so, that the form- 
factor behavior in the region of small momentum trans- 
fers is determined by the wave function at small relative 
momenta, and the behavior of the form factor in the re- 
gion of high momentum transfer — by the wave function 
at large relative momenta. The constraints for the wave 
functions provide the limitations for the form factor, and 
this is seen in the results of the calculation. 



IV. ASYMPTOTIC BEHAVIOR 

It is worth to consider the form-factor asyrnptotic be- 
havior at ^ 00 especially. In our paper [43] it was 
shown that in our approach, the pion form-factor asymp- 
totics at — > 00 and M — > does not depend on 
the choice of a wave function but is defined by the rel- 
ativistic kinematics only. We consider the fact that the 
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asymptotics obtained in our nonperturbative approach 
does coincide with that predicted by QCD as a very sig- 
nificant one. Our approach occurs to be consistent with 
the asymptotic freedom, and this feature surely distin- 
guishes it from other nonperturbative approaches. 

Let us note that it is obvious that at very high momen- 
tum transfers the quark mass decreases as it goes to zero 
at the infinity. Our approach permits to take into ac- 
count the dependence M{Q'^) beginning from the range 
where this becomes necessary to correspond to experi- 
mental data. It is possible that this will take place at the 
values of lower than 6 GeV^. 

The correct asymptotics is the consequence of the fact 
that the relativism is an intrinsic property of our ap- 
proach. To demonstrate how it works let us consider 
the simple example of point-like quarks and model wave 
fimctions (|4]). In this case in the Eqs. ([1]), ([2]): 

g^(Q') + g|(q2)^ 

For the model ^ it is easy to obtain the non- 
relativistic integral representation of the form-factor as 
the corresponding limit of the Eq. ([T|) . Now the integra- 
tion can be performed analytically and the following form 
for the nonrelativistic pion form factor can be derived: 

One can see that in the non-relativistic case, the form 
factor does not depend on the mass of constituents and 
its asymptotics can not agree with that of QCD. The 
correct asymptotics is provided by relativistic effects. 

In the relativistic case the results for the integrals can 
not be obtained analytically. To derive the asymptotic 
behavior in question it is possible to use the asymptotic 
series for double integrals obtained in Ref. ^5Q] . The first 
two terms give: 

Let us take in ([7]) the limit at M/b 0. This means 
that the parameters of the model are such that M/b ^ 1. 



The physical meaning of this limit is that the increase of 
the momentum transfer is followed by the "undressing" 
of the constituent quarks and its transformation into cur- 
rent quark of pQCD. In this limit we obtain from Eq. ([7]) 
up to logarithmic prefactors the power-like behavior co- 
inciding with that of pQCD 6] : 



V. SUMMAR 



To conclude, we make some predictions about the re- 
sults of the future JLab experiments on the pion form fac- 
tor based on the method proposed in our papers earlier. 
The method is a variant of composite quark model in the 
instant-form of Poincare invariant quantum mechanics. 
It is shown that our approach has certain advantages as 
compared with other CQM calculations. From a theoret- 
ical point of view, these advantages are the consequence 
of the fact that our approach has dispersion-relation mo- 
tivated foundations. This provides, in particular, weak 
model dependence of the results of calculations. The 
approach has demonstrated earlier its predictive power 
in describing all the data on the pion form factor ob- 
tained later in JLab experiments. Our calculations also 
give the accurate values of the pion MSR and of the de- 
cay constant /jr, and the correct asymptotic behavior at 
g2 oo. 

We hope that our model will provide a good description 
of the future JLab experiments on the measurement of 
the pion form factor in the range of momentum transfers 
up to « 6 (GeV)^. If it will occur that beginning 
from some values one needs to adjust the calculations by 
introducing the quark mass dependence on to agree 
the future data, then we propose to identify the effect 
with the appearance of perturbative effects. 

This work was supported in part by the Russian Foun- 
dation for Basic Research (Grant No. 07-02-00962). The 
work of A.K. was supported in part by the Program "Sci- 
entific and scientific-pedagogical specialists of innovative 
Russia" (Grant No. 1338). 
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